Deconvolving the information from an imperfect spherical gravitational wave antenna 
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We have studied the effects of imperfections in spherical gravitational wave antenna on our ability 
to properly interpret the data it will produce. The results of a numerical simulation are reported that 
quantitatively describe the systematic errors resulting from imperfections in various components 
of the antenna. In addition, the results of measurements on a room-temperature prototype are 
presented that verify it is possible to accurately deconvolve the data in practice. 

PACS numbers: 04.80.Nn, 95.55.Ym 

A spherical gravitational wave antenna was first proposed by Forward in 1971 The large cross section and the 
ability to determine the direction of a gravitational wave have been known for some time B, but it was not until 
recently that its potential advantages have been actively studied ^-^]. For the past few years, several groups have 
been researching the possibility of constructing large scale spherical antennas || . 

To realize the full potential of these detectors, we must first be able to understand their behavior and properly 
interpret the data they will produce. The technique needed will depend on the particular type and arrangement of 
sensors used to detect the motions of the sphere. We have proposed and investigated the "truncated icosahedral" 
(TI) arrangement |3| for six identical radial motion resonant transducers. We have shown that this arrangement on 
a perfect sphere, but with equally noisy transducers, has the property that all five quadrupole components of the 
gravitational strain tensor are measured with equal sensitivity 0. Other arrangements of transducers have been 
proposed ^|| , but most of them do not share this property JT^] . 

The TI model, also called a truncated icosahedral gravitational wave antenna (TIG A), has been criticized on the 
basis that normal (small) imperfections will break the perfect symmetry, possibly having drastic effects on such a 
degenerate system, thereby making it difficult to extract the desired information from the data. We report here the 
theoretical and experimental results of an investigation into the effects of imperfections in a spherical antenna that 
refutes these criticisms and shows that the symmetry breaking can be easily handled. 

The normal modes of an uncoupled elastic sphere can be described by the spacial eigenfunctions ^ n [ m (r, 9, (f) 
(functions of the spherical harmonics) and the purely time dependent mode amplitudes a n i m (t) ||],|Il],[l2| . In general 
relativity only the t — 2 quadrupole modes strongly interact with a gravitational wave so for the remainder of this 
discussion we will assume n = 1, £ = 2 and drop these subscripts from our notation. We denote the orientation 
of the 5 quadrupole modes relative to the lab frame by 2 Euler angles (3 m and j m . If the quadrupole modes are 
degenerate these angles are arbitrary (restricted only by orthogonality of the modes). Imperfections of the sphere can 
lift the degeneracy causing the modes to "fix" themselves in a particular orientation. This effect has been shown to be 
the main result of symmetry breaking in an uncoupled sphere; the "shape" of the quadrupole mode remains mostly 
unchanged |I^ , |l4| ]. 

We consider that a number of small resonant transducers, with index j = 1 ... J, are added to the surface of the 
sphere to mechanically transform the small motion of the surface into a large motion of the resonator. We assume 
the resonators are constructed to obey a one-dimensional harmonic oscillator equation. The particular arrangement 
of resonators is indicated by their locations {9j, 4>j) on the sphere surface and their coupling to the surface motion 
in different directions tj. In the case of ideal radial resonators we would set €j = f. The relative amount of surface 
displacement coupled to each resonator for a particular mode is specified by a "pattern matrix" defined as 
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where a depends on the material and the radius R of the sphere |2jj7|] . Again, we assume that the transducers strongly 
interact only with the quadrupole modes of the sphere, which is appropriate if the mass of a resonator is small 
compared to the mass of the sphere [||jlj| ■ 

The complete set of coupled equations of motion for the sphere and resonators can now be written [|16| 
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where we have defined 4 diagonal matrices whose elements are the mass m and spring constants k of both the sphere 
modes (superscript s) and the resonators (superscript r), 

Mj m = Sj m m m , Mj m = Sjmirij, Kj m = 5j m k m , Kj m = Sj m kj. (3) 

The column vector a has five components, one for each sphere mode amplitude, and has a one-to-one correspondence 
with the quadrupole components of the gravitational wave strain tensor ]7];[l2| . The column vector q has J components, 
one for each resonator displacement, and are the observable quantity. The column vector / contains any forces (such 
as a gravitational wave) acting on the sphere and resonators. 

A normal mode solution of these coupled equations is possible using standard techniques |l6} | . We find the response 
of the normal modes r}(u>) to any external forces, including gravitational waves, is given by 



tj(lu) = G{uj)V- 
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The matrix V is constructed from the eigenvectors of the normal modes, and the matrix G(u>) is constructed from 
their eigenvalues. The matrix V can also be used to return to the original sphere and resonator coordinates 



a(u) 
q(uj) 



= Vt,(lo). 



(5) 



V is always invertible thus it is possible to use the inverse of Eq. (||) to transform the data to normal coordinates 
where the frequency response is simple. 

Examining Eqs. and (^J), we see that in order to obtain all the information about an excitation we must have 
knowledge of both the sphere a and resonator q modes even though only the resonator motion is observable. A 
simple solution emerges if we limit ourselves to the case of "close" to the TI arrangement (six radial resonators placed 
at the center of six non-antipodal pentagon faces of an imaginary TI concentric to the sphere). The eigenfunctions 
of an ideal TIGA are such that the motion of the resonators mimic the ellipsoidal (quadrupole) deformation of the 
sphere's surface either in phase or anti-phase 0. We can picture the collective motion of the six resonators to describe 
six "resonator ellipsoids" , five of which are mimicking the "quadrupole ellipsoids" of the sphere (the sixth resonator 
ellipsoid is just a "breathing" sphere: the six resonators are moving in unison with equal amplitude and the sphere 
surface does not move at all). Each individual resonator displacement qj represents a superposition of the point 
radial deformation of the six resonator ellipsoids at a particular location. The transformation between point radial 
deformations q and ellipsoidal amplitudes g is given by the pattern matrix B: 



Bq. 



(6) 



The resonator ellipsoids g are proportional to the sphere modes a so we sometimes refer to them as "mode channels." 
We also point out that Eq. (^) can be seen as a linear coordinate transformation between the measured point radial 
displacement and a 5-dimensional abstract vector space based on the spherical harmonics that has a one-to-one 
correspondence with the five quadrupole components of the gravitational strain tensor 0,|l5|. While not unique to 
TIGA, such a simple transformation is not possible for most other arrangements of transducers. 

We are now prepared to discuss how imperfections effect our ability to deconvolve the signal from the transducers. 
One source of error is the resonator ellipsoids may be deformed so that they no longer mimic the sphere surface. 
This could lead to errors in calculating both the matrix V and the mode channels g. Another source of error is the 
uncertainty of the values of all the system parameters. Most of the parameters can be measured, either directly or 
indirectly, but we will be limited by the accuracy of that measurement. 

The transformation matrix V can be measured by applying a continuous sinusoidal force anywhere on the sphere's 
surface at the frequency of one of the normal modes. (Note that this technique requires the normal modes to be 
non- degenerate, thus it is actually preferable to have a small amount of symmetry breaking.) The frequency response 
of the resonators will be simple because they are being driven at a single frequency. From Eq. (||) we see that the 
amplitude and phase of their response make up a single column of V. By exciting each normal mode in turn, the 
complete V matrix can be measured. The only assumption made in calculating V is that the resonators are "close" 
to the TI arrangement so that an ideal pattern matrix B can be used and the sphere modes a can be replaced by the 
mode channels g. 

The quadrupole mode orientation angles f3 m and j m can be measured directly only before the resonant transducers 
are attached [lq] . This is also true for the the mass of the resonators to J and the sphere modes rn s m: and their 
respective spring constants k^ and k^. Once the system is coupled together these parameters cannot be measured 
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independently. The positions of the resonators (Oj,(f>j) can be easily measured, but precise values of their couplings 
6j are much harder to determine. 

We developed a numerical simulation of a TIGA that calculates the error associated with a measurement due to 
imperfections and asymmetries as well as our uncertainty of the precise values of the above parameters. We simulated 
a gravitational wave burst excitation of an imperfect TIGA and used the above analysis techniques to estimate the 
direction of the wave from the simulated output of the transducers. The analysis assumed parameter values of an 
ideal TIGA and a measured V matrix. The "real" (perturbed) parameters of the imperfect TIGA were calculated by 
adding uniformly distributed random numbers to the ideal values. The range of values of the random number was a 
fraction (tolerance) of the ideal value (2n for the angles) . 

The results of the simulation are shown in Fig. [l[ Plotted is the solid angle estimation error Af2 and the percentage 
error a^/h for estimating the strain amplitude for a range of tolerances. Parameters not shown had errors close or equal 
to zero. Because the tolerance is expressed as a percent of the ideal parameter values, varying the magnitude of the 
ideal parameters (for example doubling the mass of the sphere) produced identical results. Also varying the direction 
and polarization amplitudes of the gravitational wave produced similar results. From these results we conclude that 
the analysis becomes unreliable only after the tolerance of all the parameters exceeds about 3%. This is certainly an 
obtainable level of precision; one expects to set the tolerances on a real antenna to better than 1%. 

Also shown in Fig. [I] are the errors due to a finite signal-to-noise ratio, calculated in a similar fashion to the analysis 
of Zhou and Michelson [51, but using the direction finding algorithm used in our measurements Note that this 
is a random error whereas the tolerance errors are systematic. Comparing the two sources of error, we find that one 
would need a signal-to-noise ratio of about 1000 in energy before the tolerance errors become significant. While one 
might hope to observe sources at this level, the most optimistic predictions lead to considerably smaller values (lTj . 

To further evaluate the feasibility of a spherical antenna, we constructed a 0.8 m diameter, Aluminum alloy, 
room-temperature prototype TIGA Six radial resonant transducers were attached to the surface in the TI 

arrangement: placed at the center of each of the six lower pentagon faces. The response of each resonator was 
converted to an electrical signal by a piezoelectric strain gauge,whose output was recorded on a high speed data 
acquisition system. The apparatus is described in detail in Ref . p8j . 

Since no laboratory source of gravitational waves exists, we chose to apply an impulsive force to the surface of the 
prototype to test the developed techniques. To quantify the test, we checked if the location of the impulse could be 
measured from the response of the resonators. 

We began by measuring the transformation matrix V as described above. Impulsive forces were then applied 
by sending a short electrical pulse to a piezoelectric shaker attached to the surface of the TI. The data from the 
transducers were recorded, transformed to normal coordinates, fit for phase and amplitude, and finally transformed to 
mode channels. From the mode channels the location of the impulse was calculated. The algorithm used to determine 



the location of the impulse can also be used to calculate the direction of a gravitational wave [g_3 16 1 . The results of 
this analysis for the various impulse locations is shown in Fig. |^. From these results, we conclude that it was possible 
to reconstruct the location of the impulses, within the accuracy of the experiment, without major difficulties. 

In summary, from the numerical and experimental results we conclude that imperfections in a spherical antenna 
can be easily handled and that the error associated with the asymmetries are small. We find that the TI arrangement 
is fairly robust to these normal imperfections. The techniques for deconvolving the data are mostly linear algebra, 
making their implementation simple in an automated data analysis system. The measurement of the transformation 
matrix V takes into account most deviations from perfect symmetry and enables the data to be transformed to a 
space where the frequency complications can be easily handled. The algorithms were tested on the prototype TIGA 
and found to be consistent with the measured results. Finally, since all the techniques described and tested can be 
applied in situ, they are directly applicable for use on a real spherical antenna searching for gravitational waves. 

We thank W. O. Hamilton for essential advice and support. This research was supported by the National Science 
Foundation under Grant No. PHY-9311731. 
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FIG. 1. The solid angle direction estimation error and the percentage error in estimating the strain magnitude as functions 
of the tolerance on the system parameters of a TIGA and as functions of the signal-to-noise ratio for an ideal TIGA. Each solid 
line corresponds to the bottom axis and represents the results of a 200 trial simulation with the indicated parameter varied 
within the corresponding tolerances. The dashed lines corresponds to the top axis and represents the results of a simulation of 
an ideal TIGA for a range of signal-to-noise ratios. 
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